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Abstract 

It is shown that a macroscopic superposition state of radiation, 
strongly interacting with an ensemble of two- level atoms, is removed 
generating a coherent state describing a classical radiation field, when 
the thermodynamic limit is taken on the unitary evolution obtained by 
the Schrodinger equation. Decoherence appears as a dynamical effect 
in agreement with a recent proposal [M. Frasca, Phys. Lett. A 283, 
271 (2001)]. To prove that this effect is quite general, we show that 
this same behavior appears when a superposition of two Fock num- 
ber states is also considered. Higher order corrections are computed 
showing that this result tends to become exact in the thermodynamic 
limit. It appears as a genuine example of intrinsic collapse of the wave 
function. 

PACS: 42.50.Ct, 42.50.Hz, 03.65.Yz, 42.50.Lc 



1 Introduction 

Recent experimental findings in mesoscopic devices [1] seem to suggest 
that, at very low temperatures, some kind of mechanism is producing 
decoherence [2]. These same experiments seems to prove that this 
mechanism is intrinsic to the device. So, it is become demanding to 
understand a way to describe decoherence at zero temperature [3]. 
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Zero temperature decoherence in the framework of dissipative quan- 
tum systems has been devised in [4]. A general description of these 
kind of quantum systems can be found in [5]. In this paper we want to 
describe another approach to decoherence, that is, decoherence pro- 
duced by the unitary evolution in the thermodynamic limit [6]. This is 
non dissipative decoherence, a first theory of which has been given in 
[7] where time is considered as a stochastic variable. Besides, in meso- 
scopic physics, some examples of non dissipative decoherence were also 
given in [8, 9]. 

Some recent approaches rely on vacuum fluctuations to generate 
such kind of decoherence [10] and, in mesoscopic devices it appears 
that electron-electron interaction plays the dominant role. But, again 
these approaches can be leaded back to dissipative quantum systems. 

Non-dissipative decoherence can play a major role as it can make 
quantum theory self-contained without requiring arbitrary separation 
between a bath and a system. Here we demand that some systems, 
in the thermodynamics limit, develop classical behavior. When such 
systems interact with other quantum systems, decoherence indeed ap- 
pears. The evolution is anyhow unitary. 

The foundation for this kind of view of decoherence is given by 
the formal analogy between the time evolution operator of quantum 
mechanics written as exp (^—i^-^j and the density matrix of a system 
being at equilibrium exp (—(3H). In the latter case the thermodynamic 
limit has a precise meaning recovering standard thermodynamic [11]. 
Here we will show that, in the case of quantum evolution, one can 
get classical states when the system is properly prepared at the initial 
time [6]. Then, when such a system interacts with another quantum 
system, decoherence can develop. Such classical states are character- 
ized by having the system following exactly the classical equations of 
motion, for given observables, without any significant deviation due 
to quantum fluctuations, in agreement with the Ehrenfest's theorem. 

It is interesting to note that a recent experiment by Haroche and 
coworkers [12], aimed at realizing a conceptual experiment on com- 
plementarity due to Bohr, proves the appearance of classicality in the 
limit of increasing photons in a cavity. The states they have pro- 
duced are those theoretically predicted, in a pioneering work on the 
Jaynes-Cummings model in the limit of a large number of photons, by 
Gea-Banacloche [13]. In this latter work, the question of the collapse 
of the wave function is properly discussed in a similar context as ours, 
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appearing as a first hint toward the appearance of decoherence in the 
thermodynamic limit. 

In this paper we want to limit our analysis to the particular case 
of a single radiation mode interacting with iV two-level systems. So, 
by studying the thermodynamic limit, we prove that, in the strong 
coupling limit and with the two-level systems properly prepared, de- 
coherence develops destroying superposition states that are taken ini- 
tially for the field mode, driving the system toward a coherent state 
describing a classical field. Then, by computing higher order correc- 
tions, using perturbation theory, we will be able to show how these 
terms are negligible small in the thermodynamic limit, so that, the 
leading order term tends to become an exact solution in the same 
limit. This seems a clear example of dynamical collpse of the wave 
function obtained solving perturbatively the Schrodinger equation. 

The paper has the following structure. In sec 2 we introduce the 
quantum model that we want to analyze. In sec. 3 we show how the 
appearance of classical states indeed happens for an ensemble of two- 
level systems. In sec. 4 we give the leading order solution of the model 
in the strong coupling regime proving that, given a superposition of 
states, being Fock or coherent states, the system is driven toward 
a coherent state describing a classical radiation field. In sec. 5 we 
evaluate the higher order corrections to the leading order solution to 
show that, in the thermodynamic limit, our result does not change 
but rather tends to be an exact solution. Finally, in sec. 6 we give the 
conclusions with a brief discussion. 

2 The model: two-level atoms inter- 
acting with a single radiation mode 

Our aim is to discuss how decoherence can appear as a dynamical 
effect, in the limit of a very large number of systems (thermodynamic 
limit). It is a well-known matter [14] that a single two-level atom 
under the effect of an increasing number of radiation modes changes 
its behavior from coherent Rabi oscillations to spontaneous emission, 
a typical decoherent effect. 

Here, we reverse the situation by leaving a single radiation mode 
interacting with an ensemble of two-level atoms. The Hamiltonian of 
a single radiation mode interacting with a two-level atom is given by 
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[14] 



H = coa'a + — a z + go x (a + a*) 



(1) 



being u> the frequency of the radiation mode, A the separation between 
the levels of the atom, g the coupling between the radiation and the 
matter, a x and a z the Pauli matrices and a and a) the annihilation 
and creation operators of the radiation mode. In a recent experiment 
on a small Josephson-j unction circuit irradiated with microwaves [15], 
the strong coupling regime for this Hamiltonian has been obtained 
describing Rabi oscillations as theoretically described in [16]. So, the 
strong coupling regime has been practically realized in this physical 
situation. 

This model can be easily generalized to an ensemble of N two-level 
atoms as 



and again we consider the strong coupling regime as, in this case, deco- 
herence is produced dynamically by unitary evolution in the thermo- 
dynamic limit. This model has been recently studied, in a similar con- 
text, in Ref.[17]. It can be realized by an array of Josephson-j unction 
circuits, irradiated by microwaves, working as qubits for quantum 
computer. 

Let us note, at this point, that having our Hamiltonian the form 



we are neglecting K^ath^ that is, a dual situation with respect to the 
one considered by Haake and coworkers [18] that neglect H sys tem- But, 
while they describe dissipative decoherence, we aim to show, using 
perturbation theory, that decoherence is produced dynamically in the 
thermodynamic limit. 




cr z i + 9<?xi{a + aJ) 



(2) 



H — Hsystem ~\~ Hbaih 4~ H system- bath: 



(3) 
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3 Classical states by unitary evolution 
and thermodynamic limit 

In this section we want to limit our study to the Hamiltonian of the 
ensemble of two-level atoms 

N A 

#o = E ~o° zi ( 4 ) 

i=l z 

and we prove the following result [6]: 

For a proper set of initial conditions, the ensemble of two-level 
atoms evolves in time classically with respect to the variables T> x = 
J2iLi <?xi, = Eili °yi, = Ei=i &zi, in the thermodynamic limit. 

We want to emphasize that the system does not decohere but 
evolves in time washing out the quantum fluctuations in the thermo- 
dynamic limit. Besides, the system must be properly prepared as we 
cannot claim that when the system is in an eigenstate of H , behaves 
classically. 

To prove our result, we consider the following initial state 

N 

i^(o))=n(°<u)*+ 6 *it>i) ( 5 ) 
i=i 

with |aj| 2 + \bi\ 2 = 1. The set of coefficients Oj and b{ must be chosen 
in such a way that \ip(0)) is not in an eigenstate of Hq. The time 
evolution is determined through 

N 

Uo{t) = J{[e^\ [) t t (l\ + e~^\]) l ^ \\ (6) 
i=i 

giving 

N 

\m) = \{{ci i e i ^\ Di + he-^W)^. (7) 

i=l 

The average values are given by 

(E x ) = (ijj(t)\X x m)) = N(£h cos(At) + £ H sin(At)) (8) 

being i H = Ei=i(o*6i + ^b*)/N and ? H = i^Mh - a t b*)/N, 
numbers of order of unity. In the same way one has 

+N 2 (£ H cos(At)+& sin(At)) 2 
(9) 



t=i 
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so that, one gets finally that (A£ x ) 2 oc TV and the mean values are 
overwhelming large with respect to quantum fluctuations in the ther- 
modynamic limit. A similar argument runs for Ti y and T, z compo- 
nents. This means in turn that the classical equations of motion from 
the Ehrenfest theorem 

(t x ) = -A<E y > (10) 
<£„> = A(E X ) 

(£ z ) = o 

are obeyed without any significant deviation due to quantum fluctua- 
tions in the limit of very large N. This proves our assertion. 

A relevant question that can be asked with this result is if such a 
model, with this kind of initial state, can produce decoherence inter- 
acting with other quantum systems. We will answer this question in 
the next sections. 



4 Decoherence generated by unitary 
evolution: Leading order solution 

We want to show as, an ensemble of two-level atoms, strongly interact- 
ing with a single radiation mode, produces decoherence dynamically 
in the thermodynamic limit. For our aim, we start initially with a 
macroscopic quantum superposition state (generally known in litera- 
ture as a Scrhodinger cat states) given by [19] 

\^))=N{\ae^) + \ae~^))\x) (11) 

being \x) the contribution of the two-level atoms to be specified fur- 
ther, a a real number, (f> a phase and N a normalization factor given 
by 

J\f* = > (12) 

2 + 2 cos[a 2 sin(20)] exp(— 2a 2 sin 2 (f>) ' 

accounting for the nonorthogonality of the two coherent states. 

Then, we consider the strong coupling regime of the Hamiltonian 
(2) as done in Ref.[16] by assuming as an unperturbed Hamiltonian 

N 

H u = Loa)a + g^o xi (a + a)). (13) 
i=i 
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We can compute the unitary evolution operator by treating formally 
= Z)i=i °xi as a c-number obtaining 

U F (t) = e m e -i^at exp [ / 3( t ) a t _ / g(t)t a ] j (14) 

being 

|(t) = ^f-( w t-sin(c^)) (15) 



and 

P(t) = — (1 - e iw *) (16) 



operators for the ensemble of two-level atoms. 

At this stage, we choose a simplified form of the initial state (5) by 
simply taking the eigenstate |1) = ^(| |) + | |)) of a x , with eigenvalue 
1, for each atom in the ensemble, so 

ix>=nii>i- (17) 

i=i 

This kind of "ferromagnetic" state simplifies the computation giving 
us at the leading order in the strong coupling perturbation series [20] 

|V>(<)> ~ U F (t)\if;(0)) = e^WjS/'(e^W|^(t)e-^+ae i *- <wt )+e < * , W|^(t)e- i<,rt +ae-^- i<,rt ))|x). 

(18) 

We have used the property of the displacement operator for coherent 
states so to yield 

N 2 q 2 

£(*) = -^f-M " smM)), (19) 
N 9„ Jut 



P(t) = — (1 - O (20) 



and 



Mt) = -i^[P(t)e- i<t, -(3*(t)e% (21) 

<A2(t) = -if[/3(t)e^-/3*(t)e-^], (22) 

with the phases (f>i(t) and <fo(*) generated by the multiplication of 
two displacement operators. Then, it is straightforward to impose 
the thermodynamic limit N — > oo, keeping a fixed, to verify that 
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the macroscopic quantum superposition state is driven to the clas- 
sical state \(3(t)e~ tujt ), proving our assertion. The system decoheres 
removing the superposition of the states. This state describes a clas- 
sical radiation field as we are in the thermodynamic limit and so, the 
quantum fluctuations can be neglected with respect to the mean value 
[21]. 

One may ask if this is true decoherence. Actually, looking at the 
state (18) it appears as if we have just displaced the initial coherent 
state and it seems that the quantum effects are all there yet and we 
have not disposed of them as normally happens by a true decoherent 
model. We prove now, by computing the Wigner function of the state 
(18), that, in the thermodynamic limit, quantum effects are washed 
away and decoherence is recovered. Indeed, it is straightforward to 
obtain just for the interference term containing negativity and oscil- 
lations as 



Wint 



■ cxp 



7T 



x exp 



X COS 



- + ^^ (1 — cos(u;t)) — V2acos(4>) cos(u;t) 

( V2Ng . . r . \ 

p H sm{Lot) + v 2a cos((p) sin(wt) 



(23) 



2v / 2asin(0) (psin(o;t) — xcos(wt)) + a 2 sin(20) + 8a^- sin(^)(l — cos(cji)) 

id 



It is easily realized that it depends also on time that appears in 

strongly oscillating terms like cos (&a^(l — cos(wi)) sin(^)^ and sin ^8a^(l — cos(wt)) sin(0)^ 
and a sense should be attached to such terms in the thermodynamic 
limit N — > oo. Indeed, this is mathematically possible and such terms 
can be assumed to be zero in such a limit, e.g. in the sense of Abel or 
Euler [22]. A physical meaning should be attached to the oscillating 
part taken to be zero in the thermodynamic limit. Indeed, One can 
see that the time scale of variation of the oscillating part becomes 
even more smaller as ./V increases. This means that one could reach, 
in principle, oscillations on a time scale of the Planck time but one 
should expect that an average in time happens largely before this can 
happen. This is the question of the singular limits that are at the 
foundations of the theory of decoherence as pointed out by Berry [23]. 
So, the interference term can be neglected and true decoherence hap- 
pens. This point of view has also been emphasized in Ref. [8] for a 
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spin interacting with a spin bath. So, we recognize commonality be- 
tween these models and the limit in the Abel or Euler sense translates 
into an average in time. 

We want to verify if this same effect happens for a superposition 
of number Fock states. So, we consider an initial state given by 

W0)> = ^(|0) + |fc»| X > (24) 
being k an integer, yielding 

Wit)) « U F (t)\m) = e*®-j= (|/3(t)e-^> + e- ik ^\k,(3(t)e- iwt )) \x) 

(25) 

being 

oo 

|fe,/3(t)e— *) = \n)(n\D[P(t)e-^]\k) (26) 

n=0 

a displaced number state [24]. From Ref.[24] we derive the following 
formulas 

(n\D[a]\k) = [^f a^e~^Lt k \W\ 2 ) (27) 
for n > k, otherwise one has 

(n\D[a]\k) = ^y\-an (k ' n) e-^L(t n) (\a\ 2 ) (28) 

being L^x) the associated Laguerre polynomials [25]. We put 

j3'{t) = (3(t)e- iujt (29) 
and then rewrite eq.(25) as 

W)) « eW-L(\(3Xt))+ (30) 
fe-i - •- i 

— iktot 



X>> g)"[-^W] (fc -" ) e- M ^Lr«)(|/3'(t)| 2 ) + 
e -*ut g |n> ^^ (t)]( ^ fc)e _j^ L („_ fc)(|/3 , (t)|2) \ |x) _ 
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In the thermodynamic limit we can take (3'(t) to become increasingly 
large being proportional to N. So, we can approximate the associated 
Laguerre polynomials as L^(x) « for x — > oo and substituting 

into eq.(31) one gets 

(31) 

and it is not difficult to recognize the second term on the r.h.s. being 
again the coherent state \[3'{t)). Then, the radiation field is driven 
again toward a coherent state describing a classical field, the initial 
superposition state being washed out in the thermodynamic limit. 

In both cases we are left with corrections to the normalization 
factor due to the approximations we have introduced. They will not 
play any role when the average values are computed as we have applied 
a limit procedure on the wave function implying a recomputation of 
the normalization factor. 

5 Decoherence generated by unitary 
evolution: Higher order corrections 

In this section our aim will be to compute higher order corrections 
to the results in sec. 4 to prove that indeed, these corrections do not 
modifiy our argument in the given approximations. 

The dual Dyson series can be formally written as [20] 

U(t) = U F (t)Texp (-i J*dt'U F (t')H U F (t'^J (32) 

being T the time ordering operator, giving at the first order 

U {1) {t) = -iU F {t) dt'U F {t')H U F (t') (33) 

J to 

and, for the second order 

U^(t) = -U F (t) f dt'U F (t')H U F {t') f dt"U F (t")H U F (t"). 

J tg J to 

(34) 
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Our aim is to evaluate both the contributions. The first order solution 
has been put forward, firstly, in Ref.[16]. We introduce a quite general 
way to evaluate higher order corrections by introducing the reduced 
unitary evolution as 

U Rx (t) = ( X \U(t)\ X ). (35) 

We consider only the transitions with the spin bath untouched, the 
most relevant for our analysis. We will justify a posteriori this def- 
inition. This describes an unitary evolution for the radiation field. 
As we have seen in Sec. 4, the unitary evolution, at the leading order, 
produces the disappearance of superposition states in the thermody- 
namic limit. This effect, also true for the reduced unitary evolution, 
disappears for all the other orthogonal states of the spin bath that 
appear during the evolution. 

Indeed, by eq.(33) the first order gives 

(36) 

with 

-W = ^, (37) 

not dependent on N. But the first order correction modifies the spin 
bath state as 

|X') = | - l)l|l>2 ■ ■ ■ |1)JV + |l)l| - 1)2 ■ ■ ■ |1>JV + - ■ - + |l)l|l>2 ■ ■ ■ | - l)iV, 

(38) 

that is a non normalized state orthogonal to |%). We normalize this 
state by multiplying it by the factor -7= and call it \x')- The con- 
clusion is that the first order term is zero and does not contribute to 
the reduced unitary evolution. But, if we want to know the reduced 
unitary evolution with respect to the state \x') we can see that now 
is the leading order that does not contribute being zero, and we have 

uHUt) = (x'\U (1 \t)\x) = (39) 



T ^' e 4i ( JV - 1 )^( w *'- sin M')) e a(t')at-a(t')*a_ 
JO 



2 



11 



This term, in the thermodynamic limit N — > oo, has the integral wildly 
oscillating and giving a contribution proportional to jj. So, this term 
that modifies the unitary evolution as due to a modified state of the 
spin bath, can be safely neglected in the thermodynamic limit, being 
O (^7^) • This proves a posteriori, as promised, our statement on the 
relevance of the reduced unitary evolution just with respect to the 
state |x). Other non orthogonal states appear at higher orders, as we 
will see, but their contribution is even more smaller in the thermody- 
namic limit. 

We can now evaluate the second order contribution to the reduced 
unitary evolution that is given by 

uSUt) = (x\U (2 Kt)\ X ) = (40) 



2 2 

i^-§-(u)t-sin(u)t)) -Wat 8(t)a^ -/3{t)* a 



e 

N — /* dt ' e - 4i ( iV - 1 )^M'- sin M')) e -a(t')" t +a(f)*a x 

4 Jo 

r? 



f rft // e 4i(JV-l)5(art"-sin(«t")) e a(t")a t -a(*")'a 
JO 

and we have two wildly oscillating integrals. Again we can estimate 
this term, in the thermodynamic limit, as being O (^j^j an d then, it 
is even much faster in going to zero with respect to the first order 
term. This proves our assertion that, in the thermodynamic limit, the 
unitary evolution is able to wash out quantum superposition states, if 
the initial state of the spin bath is properly set. Besides, we point out 
that also the second order term of the Dyson series changes the state 
of the spin bath due to the product of Hq terms. We do not report 
this modified state here as not essential for our discussion. 

The main conclusion of this section is that, in the thermodynamic 
limit, the leading order unitary evolution tends to an exact solution 
of the Schrodinger equation. This solution, in the same limit, can 
produce decoherence washing out quantum superposition states. 



6 Discussion and Conclusions 

We have showed that a single radiation mode, interacting with an 
ensemble of two-level systems, can undergo decoherence during its 
unitary evolution. The proof is obtained using perturbation theory 
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and assuming strong coupling between the field and the ensemble of 
two-level systems. At the leading order, the unitary evolution washes 
out superposition states when the thermodynamic limit is taken and 
the ensemble of two-level systems is properly initially prepared. The 
proof is given both for a macroscopic quantum superposition state and 
for a superposition of Fock states. Higher order corrections till sec- 
ond order are then computed, proving that the leading order solution 
tends to be an exact one in the thermodynamic limit. So, the classical 
radiation field is an attracting solution when a radiation mode inter- 
acts with a very large ensemble of two-level atoms. It is interesting to 
note that this solution of our model seems to be a genuine example of 
dynamical collapse of the wave function. 

This result opens up the possibility to prove a quite general result. 
That is, when N non interacting particles can realize classical evolu- 
tion of observables, one may ask if, interacting with such a system and 
assuming the proper initial state, such an effect of decoherence may 
be an ubiquitous one. Indeed, our approach proves to be a useful ex- 
tension of the formal analogy between statistical physics and quantum 
mechanics. 
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